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ABSTRACT 
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Introduction 


In the past few years there has been a rapid development of computational 
methods for identification and control of systems governed by functional and 
partial differential equations. The recent literature on this topic is 
extensive and we will not attempt a review of the area. The purpose of this 
paper is to show, via a fairly simple example, that care must be taken to 
ensure that finite dimensional approximations of distributed parameter systems 
preserve important system properties (i.e., controllability, observability, 
stabilizability, detectability, etc.). Moreover, if the particular scheme 
used to construct the finite dimensional model does not take into account 
these system properties, then the model may not be suitable for control design 
and analysis. 

Clearly, controllability and observability properties of the finite 
dimensional models depend on both the distributed parameter system and the 
type of approximation scheme used to construct the lumped model. It is also 
important to note that although a particular numerical scheme may be 
convergent (i.e., consistent and stable) and well suited for simulation of 
open-loop dynamics, the scheme may not be suitable for use in the design and 
analysis of feedback control systems. A finite dimensional model to be used 
in a LQR type design algorithm should be based on a numerical scheme that 
approximates both the distributed parameter system and its adjoint (see [2], 
[5], [7], [8]). The basic ideas can best be illustrated by a simple example. 

A Hybrid Control System 

Hybrid control systems are systems governed by coupled partial and 
ordinary differential equations. Such systems occur in large flexible 
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structures and are typical models for boundary control problems when actuator 
dynamics and/or tip masses are included in the analysis. We shall focus on 
the simple cable-spring-mass system shown in Figure 1. The cable of length L 
has lineal density a and is under tension t. The device at the right end 
of the cable maintains the tension and provides no impedence to the vertical 
motion of the mass m. Let y(t,x) denote the displacement of the cable from 
the x— axis and u denote an applied force acting on the mass. 

The energy of the mechanical system (Kinetic plus potential) is given by 

E(t) = / {ay 2 (t,x) + xy 2 (t,x)}dx + ky 2 (t,L) + my 2 (t,L). (1) 

0 c 

for 0 t < + 00 . There are several possible state space models of this 
system. We shall consider the approach followed in [3] and formulate the 
problem in the state space 

z = L^COjL) x L2(0,L) x It x 13, (2) 

If Zj = col(<J> j ,<}>2 ,n i ,ri 2 ) and Z 2 = col(i|)^ ,^ 2 ,5 ^ >^ 2 ^ > then the inner product 
defined by 

L L 

< Z 1 » z 2> = / T* 1 (x)* 1 (x)dx + / a<f> 2 (x)ij)2(x)dx + kn^ + m 2 K 2 (3) 

leads to the energy norm |z|3 = <z,z> for the system described in Figure 
1. Moreover, if we identify z^t.x) = y x (t,x), z 2 (t,x) = y t (t,x), 

= y(t,L) and = p^(t) = y^(t,L), then the system can be 

described by the hybrid system (here a? = t/o > 0): 
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with boundary conditions z 2 (t,0) = 0 and z 2 (t,L) = p 2 (t). 
Define the operator A on the domain 


D(A) = {z - coK+j, ♦ 2 ,n 1 ,n 2 ) |* ± e W 1,2 (0,L), i-1,2 

and <f> 2 (0) = 0, <t» 2 (L) = n 2 } (5) 
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where "prime" denotes The operator B:R-»-z is defined by 


( 7 ) 


With A and B as above the control system in Figure 1 may be realized as 


2(t) = Az(t) + Bu(t) (8) 


with initial data 



z(0) = z Q = col(<(. 1> <j) 2 ,n 1 ,Ti 2 ). (9) 

It can be shown that A generates a Cp-semigroup on Z and for suitably 
restricted initial data solutions to (8) — (9) are equivalent to solutions of 
the system 


y tt (t,x) = a 2 y xx (t,x) 
my tt (t,L) = -ky(t,L) - ty^t.L) + u(t). 


( 10 ) 


Note that (formally) the system (4) is the symmetric hyperbolic form of (10). 
In order to define an output operator we let 


. L/k + e 

= 2i / < l ,(x)dx 

L/k - e 
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denote the mean-value operator at L/k. Define 


C :Z+tH by 
e 


♦i 

♦2 

n 2 


CjMjjC + j) 

C 2 V*2> 

C 2 M 2 ( *2> 


C 2 n 2 


C 3 n l 


(ID 


-3 

where e > 0 is small (e < 10 is sufficient for the numerical results 

presented below) and c^, i = 1,2,3 are positive constants. The operator 

acting on z(t) "observes" the "average" slopes and velocities at 

x = *1 h. and L and the displacement at x « L. 

4 5 2 

Consider the problem of minimizing 

T 

J (u) = / {lie z ( t ) II 2 + R|u(t) ) 2 }dt, (12) 

0 e 

where z(t) is the mild solution to (8) - (9) and 0 < T +°° is fixed. 
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In order to solve this optimal control problem, one must introduce some 
type of approximation scheme. We shall discuss two schemes to compare their 
use in the solution of the above optimal control problem. 

The first scheme is detailed in [3] and was used to estimate parameters 
in the system described in Figure 1. We shall give a brief description of the 
scheme. Divide the interval [0,L] into N = 2^ subintervals with nodes 
= ih/N, i = 0,1,2,...,N. Let h^(x) denote the unique continuous, 

piecewise-linear function satisfying h^(x^) = 6 . These functions are the 

so-called hat or Chapeau functions. For i = 0,1,2,...,N define u^ e Z 

by u ; = col(h^(» ), 0,0,0), and for i = 1,2, ...,N define e Z by 

= col(0,h^(» ),0,h^(L)). Let Vq denote the vector 

N i\t 

v Q = col(0, 0,1,0). Consider the finite dimensional subspace Z c Z defined 
by 

Z N = span{u^,v^, i = 0,1,2,... ,N} (13) 

and let P N denote the orthogonal projection from Z onto Z N . Note that 
the set {u^, i = 0,1,. ...N} {v*jj, i = 0,1,... ,N} is a basis for Z N 
and dim Z N = 2(N + 1). 

Moreover, it is important to note that D(A)CZ N for each 

N=2^ (k=l ,2, . . .) and hence we may define an approximating operator by 

A n = P n AP n (14) 

and approximating operators B N :R-*-Z N by B N = P N B and C N :Z N ->-l7 by 
N N N 

C z = C^Z , respectively. 
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Using these definitions one can show that the sequence of finite 
dimensional models 


z N (t) = A N z N (t) + B N u(t) 

z N (0) = zjj = P N Z Q (15) 

y N (t) - c N z N (t) 

provides a stable and consistent approximation to the distributed parameter 
system (8) - (9) with output defined by (11) • Moreover, since this problem is 
skew-adjoint, this scheme provides a stable and consistent approximation to 
the adjoint system. 

This particular scheme was used in [3] to estimate various system 
parameters in the system model. It worked quite well when used for such 
parameter estimation problems. However, when one attempts to use system (15) 
to compute suboptiraal gain operators by minimizing the functional 

T 

J N (u) = / {IIC N z N (t)ll 2 + R|u(t)| 2 }dt, (16) 

0 

then a number of problems occur (both theoretical and numerical). First it is 
easy to see that the symmetric hyperbolic/hybrid equation (4) with boundary 
conditions 


z 2 (t,0) = y t (t,0) = 0 and z 2 (t,L) = y t (t,L) = y 2 (t) 


does not take into account the (physical) constraint y(t,0) = 0. In 

particular, there is a "rigid translation" mode present in the distributed 
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parameter model (8) - (9) that leads to a lack of controllability in the 
finite dimensional model (15). However, there are other problems with this 
particular scheme that lead to the loss of controllability/observability 
properties and effect the performance of numerical algorithms based on the 
model (15). Note that the symmet ric-hyberbolic form (4) of the wave equation 
(10) requires that a condition of continuity be satisfied. In particular, the 
condition 


y tx (t,x) = y xt (t,x) (17) 

that is implied by (4) is not reflected in the approximating system (15). 
There should be a continuity condition relating the "elements" u^ that 

approximate y x (t,x) and the "elements" vj? that approximate y t (t,x). 
The use of finite elements that satisfy such continuity conditions has been 
proposed for various plate and beam models in order to prevent so called 
"locking" elements (see [6]). 

In order to see how these considerations can lead to "better" finite 
dimensional models for control, we construct a second approximation scheme. 
Let f^ = col ( 0 , h^ ( • ) , 0 , h^ (L ) ) and e^ - col(-^ h^( • ) ,0 ,h^(L) ,0) for 

i = 1,2,...,N. Define W N cZ to be the finite dimensional space 

W N = span{e^,f*J, i - 1,2, ...,N} (18) 

and note that dim W N = 2N. Let Q N :Z*W N denote the orthogonal projection 
onto W^. Since does not contain D(A) it is not possible to define 

by A^ = Q^AQ^. A rigorous deriviation of the approximating operator A^ 
relies on the theory of sesquilinear forms and is similar to the approach used 
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in [1] and [8], The basic idea is to expand the state z(t,x) in the basis 
elements e^ and and use the weak form of (8) - (9) to construct the 
finite dimensional model. In particular, let 


z N (t,x) 


r z^(t,x) 


N,,. \ 

z 2 (t,x) 


N,. v 


L ^(t) 


l y*?(t)e^(x) + w^(t)f^(x) 
i=l 


(19) 


and substitute z N (t,x) into the weak form of (4) (see [1], [8]) to obtain a 
finite dimensional model for the coefficients y|?(t) and w^(t). ^ 0ne 
defines x N (t) e l£ N by 


y 

y 


N 

1 

N 

2 


(t) 

(t) 


N,„* - 

x (t) = 


N,.* 

y N (t) 


Wj(t) 

w 2 (t) 


w S ( t) 


( 20 ) 
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then the resulting finite dimensional model becomes the 2N-dimensional system 
of the form 


M N x N (t) = F N x N (t) + G N u(t) 

( 21 ) 

x (0) = x n 


with output 


y N (t) = D N x N (t) , (22) 

where D N is a matrix representation of the operator C restricted to 

£ 

W N , and Xq is the vector of coefficients of the (representation of the) 
projection of z Q onto W N . 

This particular scheme has a number of nice features that make it 
suitable for control design. First observe that in (19) one can show that 

z j(t,x) = ~ z!J(t,x) (23) 

and 

yj(t) = wj(t), i = 1,2,... ,N. (24) 

Thus, the continuity condition (17) is preserved. Moreover, the expansion 

(19) removes the uncontrollable "rigid translation" mode that appears in the 
model (15) and numerical checks show that the finite dimensional model (21) is 


controllable. 
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Control synthesis was attempted with both finite dimensional models. The 
numerical algorithm used for computing the (time-varying) feedback gain 
operator is based on Chandrasekhar factorization of the Riccati operator 
equation (see [4], [9]). The system parameters and initial conditions are the 
same as in [3]. The time interval is [0,48 sec], the output weights are 
C^ = 50., C 2 = 25., C 3 = 1 . and the control weight is unity. 

Several numerical experiments were attempted with the design model 
(15). Various values of the grid parameter N were tried, system parameters 
were varied and we even added viscous damping to the model. None of the 
results were satisfactory. 

On the other hand, numerical results with the new model, described by 
equations (19) - (22), were quite good. Shown in Figures 2-5 are 
comparisons of open-loop and closed-loop responses for both y x (t,x) and 

y t (t,x) at several locations along the cable. These results were obtained 
with grid parameter N = 16. Several simulations were run with the gain 
operator fixed at its steady-state value. To the scale shown, the graphs were 
not distinguishable from those presented here. 

Note that both schemes are convergent to the original distributed 
parameter model and both finite dimensional models are sufficient for various 
simulation and identification studies. However, the important point is that 
the second scheme preserves a number of system properties that make it more 
suitable for control design. 
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Figure 2: Time History of v 16 (t -) 
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Figure 5: Time History of y^(t, L) 
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